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AND FLUCTUATIONS OF RANDOM MATRICES: 
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Abstract. We extend the relation between random matrices and 
free probability theory from the level of expectations to the level of 
fluctuations. We introduce the concept of "second order freeness" 
and interpret the global fluctuations of Gaussian and Wishart ran- 
dom matrices by a general limit theorem for second order freeness. 
By introducing cyclic Fock space, we also give an operator alge- 
braic model for the fluctuations of our random matrices in terms of 
the usual creation, annihilation, and preservation operators. We 
show that orthogonal families of Gaussian and Wishart random 
matrices are asymptotically free of second order. 



1. Introduction 

Free probability has at least three basic facets: operator algebras, 
random matrices, and the combinatorics of non-crossing diagrams. This 
can be seen very clearly in Voiculescu's generalization of Wigner's semi- 
circle law to the case of several independent matrices jVoilj . The dis- 
tribution arising in this limit of random matrices can be modelled by 
a sum of creation and annihilation operators on full Fock spaces and 
described nicely in terms of non-crossing partitions. 

On the random matrix side, Voiculescu's theorem describes the lead- 
ing contribution to the large iV-limit of expectations of traces of Gauss- 
ian random matrices. However, in the random matrix literature there 
are many investigations on more refined questions in this context. On 
one side, subleading contributions to the large iV-limit are of interest 
and have to be understood up to some point for dealing with ques- 
tions concerning the largest eigenvalue of such random matrices. On 
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the other side, there has also been a lot of interest in leading con- 
tributions to other important quantities, like, e.g., global fluctuations 
(i.e., variance of two traces) of the considered random matrices. One 
should note that there are relations between these two questions. We 
are not going to explore these relations here, but we want to direct 
the reader's attention to the so-called "loop equations" in the physical 
literature (see, e.g. |Eyn| ) and to the "master equation" in |HTj . 

We will concentrate in this paper on the second kind of question. As 
is well-known from the physical literature, in many cases these leading 
contributions are given by planar (or genus zero) diagrams and thus 
have quite a bit the flavour of the combinatorics of free probability. In 
the recent paper |MNj this description was made precise for the global 
fluctuations in the case of Wishart matrices, and in particular the rele- 
vant set of planar diagrams ( "annular non-crossing permutations" ) was 
introduced and examined. However, this description of the fluctuations 
in the large iV-limit was on a purely combinatorial level. Since it looks 
quite similar to the description of free Poisson distributions in terms 
of non-crossing partitions, one expects to find some genuine free prob- 
ability behind these results. In particular, one would expect to have a 
description on the level of operator algebras and to have also a precise 
statement of the kind of "freeness" that arises here. 

In this paper we will show that this is indeed the case. On one hand, 
using the notion of a cyclic Fock space, we can formulate the fluctu- 
ations in terms of the usual creation, annihilation, and preservation 
operators. On the other hand, we will also introduce an abstract "free- 
ness" property for bilinear tracial functionals, which not only give us 
a conceptual understanding, but, on the other hand, is also crucial for 
proving our main theorems on the fluctuations. 

Second order freeness, while stronger than the freeness of Voiculescu, 
nevertheless appears to be a central feature of ensembles of random 
matrices. Indeed, in this paper we prove that two standard exam- 
ples of random matrix ensembles exhibit second order freeness: or- 
thogonal families of Gaussian random matrices and orthogonal families 
of Wishart random matrices are asymptotically free of second order. 
Moreover in MSB] we show that independent Haar distributed random 
unitary matrices are asymptotically free of second order. 

The main results of the paper are thus. In section 5 (with proofs 
in section 7) we show that semi-circular and compound Poisson fam- 
ilies on the full Fock space are free of second order. In section 6 we 
establish the basic properties of second order freeness and prove a gen- 
eral limit theorem. In section 8 we diagonalize the fluctuations in the 
Gaussian and Wishart case, thus recovering and extending results of 
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Cabanal-Duvillard |C-Dj . In section 9 we prove asymptotic freeness of 
second order for orthogonal families of Gaussian and Wishart random 
matrices. 

2. Preliminaries 

Here we collect some general notation and concepts which we will 
use in the following. 

Our presentation should be, by and large, self-contained, however, it 
will rely of course on the basic ideas and concepts of free probability. 
For more details on this, one should consult |VL)JN|. I Voi2|. |NSp[ HP . 
Furthermore, the concepts of annular non-crossing permutations and 
partitions will play a crucial role. We will provide all relevant infor- 
mation on them in the text. However, our presentation will be quite 
condensed, and for further details one should consult the original paper 
jMNj. 

2.1. Some general notation. For natural numbers n, m G N with 
n < m, we denote by [n, m] the interval of natural numbers between n 
and m, i.e., 

[n, m\ := {n, n + 1, n + 2, . . . , m — 1, m}. 

For a matrix A = (ay)fy =1 , we denote by Tr the un-normalized trace 
and by tr the normalized trace, 

- 1 
Tr(A):=5>,, tr(A) := —Tr(A). 

i=i 

For anneN, we will denote by P(n) the set of partitions of [1, n], 
i.e., a = {B\, . . . , B r } G P{n) is a decomposition of [l,n] into disjoint 
subsets Bf. Bi ^ for i = 1, . . . , r, fl Bj = for i j and 

r 

[l,n] = \jB l . 

i=l 

The elements Bi of o will be addressed as blocks of a. 

Given a mapping i : [l,n] — > [1,N], the kernel, ker(i), is defined as 
the partition of [1, n], such that two numbers k,l G [1, n] belong to the 
same block if and only if i(k) = 

If we are considering classical random variables on some probability 
space, then we denote by E the expectation with respect to the cor- 
responding probability measure and by k r the corresponding classical 
cumulants (as multi-linear functionals in r arguments); in particular, 

ki{a} = E{a} and k2{ai, = E{aici2} — E{ai}E{a2}- 
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2.2. Annular non-crossing permutations and partitions. The 

leading asymptotics of various random matrix quantities can be de- 
scribed in terms of special "planar" objects (see, e.g., |Eyn[ IZvoj ). 
There are two equivalent ways of formulating these results: a geo- 
metric "genus" -expansion, expressed by a sum over surfaces where the 
planar part corresponds then to sums over surfaces of genus zero; an 
algebraic description, where instead of using surfaces one can sum 
over permutations and planarity is then a geodesic-like condition on 
a length-function of these permutations. If one prefers to associate 
partition like pictures with permutations, then planarity is a condition 
that these partitions have non-crossing diagrams (where, however, one 
has to be careful about which drawings are allowed). 

We prefer to think in terms of permutations and partitions. Let us 
recall the relevant definitions and results. 

Let, for r > 1, natural numbers n(l), . . . ,n(r) be fixed. Consider a 
partition a G P(n(l) + - ■ - + n(r)). In |MNj . the class of "multi-annular 
non-crossing partitions" NC(n(l), . . . , n(r)) was defined and, for r = 2 
("annular" case), an extensive study of various characterizations of 
such non-crossing partitions was made. We will not go into details 
here, but we only want to state the characterization which we will use. 
It will be the case r = 2 which is relevant for us; so let us use the 
notation NC (n, m) for the (n, m)-annular non-crossing partitions. It 
is a good picture to think of two concentric circles, with n points on 
the outer and with m points on the inner. We put [l,n] in clockwise 
order on the outer circle and [n + 1, n + m] in counter-clockwise order 
on the inner one. Adopting this convention will require that in some 
of our formulas the indices corresponding to the outer circle run in the 
opposite direction from the indices on the inner circle. 
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Figure 1. On the left is the annular non-crossing permuta- 
tion (1, 2, 12, 9, 8) (3, 4) (5, 10, 11) (6) (7). On the right is the 
permutation after the blocks that are contained in one circle 
have been removed. 

Consider a a G P(n + to). We shall give a recursive procedure for 
deciding if a is annular non-crossing. Suppose a has a block which 
is contained either in [l,n] or in [n + l,n + to] and which consists of 
cyclically consecutive numbers; then we remove this block and repeat 
the process until we get a partition a' G P(n' + to') with no blocks 
which are contained in either [1, n] or [n + 1, n + to] and which consist 
of cyclically consecutive elements. Then, by definition, a will be in 
NC(n,m) if and only if n',m' > 1 and a' G NC(n' ,m'). Thus it 
suffices to say when a is in NC (n, m) for a with no blocks which are 
contained in either [1, n] or [n+1, n+m] and which consist of cyclically 
consecutive elements. 

The characterizing property of such a is the following: If we write 
the blocks B G cr in the form B = B' U B" , where 5' C and 
B" C [n + 1, n + to], then, for all blocks -B of a, both parts, £?' and 
B" must be non-empty and each of them must consist of cyclically 
consecutive numbers. Furthermore, the cyclic order of the restrictions, 
B[, B' 2 , . . . , B' k , of the blocks of o to the interval [l,n] must be the 
reverse of the cyclic order of the restrictions, B'[ , B%, . . . B'J. , of the 
blocks to the interval [n + l,n + m]. Note that this characterization 
contains the statement that a a G iVC (n, m) is connected in the sense 
that at least one block of a contains elements both from [1, n] and from 
[n + 1, n + to] (i.e., a connects the two circles). 

In the context of random matrices, it is permutations, not partitions, 
which appear in calculations. In order to go over from partitions a to 
permutations it one has to choose a cyclic order on each of the blocks 
of a. Choosing such an order for each block will produce an "annular 
non-crossing permutation" out of an annular non-crossing partition. 
The set of annular non-crossing permutations is denoted by S^d 71 , Tn) 
- and by SVcWl)) • • • , n{r)) in the multi-annular case - and it is this 
set which was the main object of interest in [MN]. In [MN, §6] it 
was shown that a permutation r is in SNc( n i m ) if two conditions are 
satisfied. The first condition is connectedness: at least one cycle of 
t connects the two circles and the second is planarity: the geodesic 
condition must be satisfied: #(7r) + #(7r _1 7) = to + n, where #(vr) 
denotes the number of cycles of 7r and 7 is the permutation with two 
cycles: 7 = 7„, m = (1, 2, 3 . . . , n)(n + I,..., n + m) 
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We wish to describe what it means to be a a non-crossing permuta- 
tion on an r multi-annulus. By an r multi-annulus we mean a collection 
of r circles with n(l) points on the first circle, n(2) points on the sec- 
ond circle, . . . , n(r) points on the r th circle. Connectivity of r means 
that every pair of circles is connected by at least one cycle of r. The 
planarity of r is defined using a geodesic condition. Let 7 n (i),n(2),...n(r) 
be the permutation of [n(l) + • • • + n(r)} with r cycles — the r cy- 
cle being (n(l) + ■ • • n{r — 1) + 1, ... , n(l) + • ■ ■ n{r)). r will be pla- 
nar if t satisfies the geodesic condition #(r) + #(7n(i),n(2),...n(r) r_1 ) = 
n(l) + • • • + n(r) + 2 — r. 

As observed in jMNj . there is not necessarily a unique choice of 
a cyclic order on a block of cr; to put it another way, the mapping 
from 7i to a (which consists in forgetting the order on the cycles) is 
not injective. However, this deviation from injectivity is not too bad. 
Let us consider a block B G a, and denote by B' := B R [l,n] and 
5" := Bn [n + 1 , n + m] the parts of B lying on the first and on the 
second circle, respectively. On each of the two circles we respect the 
given cyclic order on (1, . . . , n) and on {n + 1, . . . , n + m) The allowed 
orders on B thus consist of choosing a 'first' element of B' and a 'first' 
element of B"\ then the order on B is obtained by running through B' 
from the first to the last element, then going over to the first element in 
B" and continuing in B" to the last element. Hence the only freedom 
we have is the choice of first elements in B' and in B". 

Let us call a block B G a a through- block, if both B' :— B H [l,n] 
and 5" := 5n[n + l,ri + m] are non-empty. Then only a with exactly 
one through-block have two or more 7r's in SVc( n ? m ) as preimages. 
Namely, if B = B' U B" is the unique through-block of such a <r, then 
every element from B' can be chosen as first element, and the same for 
B", thus there are exactly \B'\ ■ \B"\ possible choices of cyclic orders 
for B. If, however, there is more than one through-block, then the first 
element on each component of them is uniquely determined and there 
is exactly one possible order for each block. 

For Gaussian random matrices only non-crossing pairings will play 
a role. These are those annular non-crossing partitions for which 
each block consists of exactly two elements. One should note that 
in this case the distinction between permutations and partitions van- 
ishes, because for pairings there is always exactly one possibility for 
putting an order on blocks. We will denote the set of annular non- 
crossing pairings by iVC^n, m); and, for the multi-annular situation, 
by NC2(n(l), . . . , n(r)). In the multi-annular case the geodesic condi- 
tion can be written #(7m,...,n r 7r ) = 2 — r + (ni + ■ • • + n r )/2. 
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3. Combinatorial description of global fluctuations 

We are interested in the fluctuations of various types of N x N ran- 
dom matrices around their large N- limit. Here, we are going to consider 
two classes of random matrices, namely Gaussian random matrices and 
(a generalization of) Wishart matrices. Let us fix the notation for our 
investigations. 

3.1. Semi-circular case. Let (Xn)n&n be a sequence of N x N Her- 
mitian Gaussian random matrices. Then, in the limit iV — ■> oo, 
converges to a semi-circular variable s. Let us consider directly the 
case of several such Gaussian random matrices. The entries of different 
random matrices need not be independent from each other, but they 
have to form a Gaussian family. A convenient way to describe such a 
situation is to index the matrices by elements from some real Hilbert 
space TCm, such that the covariance between entries from Xjy(f) and 
X^(g) is given by the inner product (/, g). More precisely we say that 
{X n (f)}f e -n m is a family of Hermitian Gaussian random matrices if 
X N {f) = (x<j(/))£. =1 and the entries {x itj {f) | 1 < i,j < N,f G H R } 
form a Gaussian family with covariance given by 

¥>{xij(f)x k i(g)} = for i < j, k < I, and f,g eH R and 

E{xij(f)x k i(g)} = S ik 5ji ■ jy(f,g) for i<j,k<l, and f,g e H R 
By Wick's formula (see e.g. PJ Sections 1.3 and 1.4]) we have 

X %2 ,32 ' X i2k ,32k (/2k)} = Yl II ( X ir,3r(fr),Xi s ,j s (fs)) 

ir£V2(2k) (r,s)en 

where the sum is over all pairings it of [2k] and the contribution of each 
pairing is the product of (x ir j r , x is j a ) over all pairs (r, s) of ir. 

Then, Voiculescu's extension of Wigner's theorem to this multi- 
dimensional case states that, for iV — > oo, such a family of random 
matrices converges to a semi-circular system with the same covariance. 
We want to look more closely on that convergence and investigate the 
"global fluctuations" around this semi-circular limit; this means, we 
want to understand the asymptotic behaviour of traces of products of 
our random matrices. It turns out that, with the right scaling with 
N, these random variables converge to a Gaussian family and thus the 
main information about them is contained in their covariance. 

If one invokes the usual genus expansions for expectations of Gauss- 
ian random matrices then one gets quite easily the following theorem. 
It turns out that the leading orders are given by planar pairings. Since 
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we are looking on cumulants and not just moments, the relevant pair- 
ings also have to connect the r circles. 

Theorem 3.1. Let Xjsr(f) (f G Hr) be a family of Hermitian Gaussian 
random matrices. Let k r denote the r th classical cumulant (considered 
as multi-linear mapping of r arguments) , then for f\, . . . , f ni+ ... +rir G 
Hr, the leading order of the cumulants of the random variables 

(1) {Tr{X N {f l )---X N (f n{l) )),..., 

7V(^jv(/n(l)+-"+n(r-l)+l) ' ' • ^iv(/n(l)+-+n(r))) } 

are given by 



(2) kr\Tr[X N {f x )---X N {f n{1) )l--- , 

Tr[X N (f n {l)+-+n(r-l)+l) ' ' ' X N (f n (i) + ... +n (r))] | 

= N2 ~ r - e n (fi,fi)+o{N-*). 

n£NC2(n(l) ,.. . ,n( r )) (*j)G7r 

Proof: Let n = n(l) + ••• + n(r), 7 be the permutation of \n% + 
• • • + n r ] with the r cycles (1, . . . , ni)(ni + 1, . . . ,ni + 7^2) • • • (ni + 

h n r _i + 1, h ni H h n r ), and F; = tr(X A r(/ ri(1)+ ... +n( ^i )+ i) 

■ ■ ■X N (f n{ i) + ... +n ( i) )). Then by Wick's formula 

JV 

E(Y 1 ---Y r ) = (A) ' ' ' x wi (/n)) 

ii,...,i n =l 

iV 

»l,...,i»=l 7rg-p 2 (n) (k,l)e-rr 

N 

7re"P2(n) (k,l)en ii,...,i„=l 

= e n (hJi)N* M 

Following the argument in [MN, proof of Proposition 9.3] we have 



TreT^n) (fc,Z)€7r 
7r is connected 



The terms of highest order are the planar ones thus we obtain equa- 
tion ©. □ 
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This theorem contains all relevant combinatorial information about 
the asymptotic behaviour of our traces. Since an increase of the num- 
ber of arguments of the cumulants corresponds to a decrease in the 
exponent of N, a cumulant k r will always dominate a cumulant k p if 
r < p. So in leading order only the first cumulant survives in the limit, 
which gives us the following statement analogous to the law of large 
numbers. 

Corollary 3.2. For each fi, . . . , f n e Hm, the random variables 
{tr(x N (f 1 )---X N (fn))} N 

converge in distribution to the constant random variables a(fi, . . . , f n ) ■ 
1, where 

a(/ 1> ...,/„)= Yl II </«»/;>■ 

tt€N C 2 (n) (i,j)EiT 

This corollary is of course just a reformulation of Voiculescu's result 
that E{tr[Xjv(/i) • • -X N (f n )]} has, in the limit N — > oo, to agree with 
the corresponding moments of a semi-circular family. 

But we can now go a step further. If we subtract the mean of the 
random variables, then the first cumulants are shifted to zero and it 
will be the second cumulants which survive - after the right rescaling. 
Since higher cumulants vanish compared to the second ones, we get 
Gaussian variables in the limit. 

Corollary 3.3. Consider the (magnified) fluctuations around the limit 
expectation, 

F N (h, . . . , f n ) := N ■ {tr[X N {h) ■ ■ -X N (f n )} - a(f u . . . , /„)) 
= Tr[X N U\) ■ ■■X N (f n )\ - Naifr, fn)- 

The family of all fluctuations (F/vf/i, • • • , f n )) cM f c ^ converges in 
distribution towards (F(fi, . . . , f n )) , a centered Gaussian fam- 

ily with covariance given by 

E{F{h, ...,/„)• F(f n+1 , f n+m )} = II ^' fi)' 

weNC2(n,m) (i,j)67T 

Our goal now is to present a conceptual understanding of this form 
of the covariance; in particular one that would easily diagonalize it. 
In principle, this is a purely combinatorial problem. However, our 
point of view is that limits of random matrices which have the flavour 
of free combinatorics should also have a description in terms of the 
operator side of free probability, i.e., operators on full Fock spaces. 



10 



J. A. MINGO AND R. SPEICHER 



We will provide such a description and show that it diagonalizes our 
covariance. 

3.2. Compound Poisson case. Let (Xjv) Nm be a sequence of N x iV 
complex Gaussian random matrices (i.e. the entries of are indepen- 
dent centered complex Gaussians with variance 1/N) and let (-DTv)neN 
be a sequence of N x N non-random matrices for which a limit distri- 
bution exists as iV — > oo. Then, in the limit iV — > oo, {X^, X^, D^} 
converges in distribution to {c,c*,d}, where c is a circular element, d 
has the limit distribution of the Dm, and {c, c*} and d are free. In par- 
ticular, X n DnXn converges to c*dc, which is a free compound Poisson 
element, see |Sp2| 4.4]. We shall discuss the fluctuations of the random 
matrices 

Pn '■= X* N D N X N 

around the limit c*dc. Since P/v is a generalization of a Wishart matrix, 
we will call it in the following a compound Wishart matrix. 

As we shall see it is appropriate to consider a more general situation. 
Namely, consider not just a single non-random matrix Dn, but also all 
its powers D k N at the same time, or more generally, let us consider 
a family {D[ N \ D% N \ D^ N \ . . . , Dp N ^} N of N x N complex matrices. 
We shall say the family converges in distribution if there are operators 
di, d 2 , d 3 , . . . , dp and a tracial state if) on V, the complex *-algebra 
generated by {1, d\, d 2 , ■ ■ ■ , d p }, such that 

hm tr[D^...D^]=4j(d H ---d lk ) 

for all ii,i 2 , ■ ■ ■ ,ik- 

We are again interested in global fluctuations of these matrices in the 
limit iV — > oo; i.e., we want to consider the asymptotic behaviour of 
mixed moments of our random matrices. Again the key point is the un- 
derstanding of the leading order of the cumulants in these traces. This 
leading order is described by summing over non-crossing permutations, 
but in contrast to the semi-circular case, all permutations contribute, 
not just pairings. In order to describe the contribution of such a general 
non-crossing permutation, we need the following notation. 

Notation 3.4. Let (A,ip) be a unital algebra with a tracial state ip\ 
for each ir G S p we shall define a p-linear functional, ip n , on A x • • • x A. 

Write ir = ci • c 2 Ck as a product of disjoint cycles, and for each i, 

c i — ( r i,i, ■ ■ ■ , r i,h)- Then define the p-linear functional ip n by 

k 

^(01,02,03, ...,a p ) = Y[ip{a n>1 ■ ■•a rM .) 

i=i 
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Note that we need if) to be a trace, because a cycle c comes only with 
a cyclic order. 

An example of this notation is the following, take 

tt = {(1,2, 6), (3, 4, 5)} 6 S NC (S, 3). 

Then 

ip w (a 1: a 2 , a 3 , a 4 , a 5 , a 6 ) = ^(aia 2 a 6 ) • ^(a 3 a4a 5 ). 
Note also that the cyclic order is important. In Sjvc(2, 1) consider 
7r 1 = {(1,2,3)} and vr 2 = {(1, 3, 2)}. 

Although their block structure is the same, as permutations they are 
different elements from Snc{2, 1) an d we have 

ip 7T1 (a 1 ,a 2 ,a 3 )=ip(a 1 a 2 a 3 ) and V-V 2 ( a i, °2, a 3 ) = ip^a^). 

We shall denote the number of cycles in the permutation tt by #(7r). 

Let us now state the basic combinatorial description of the leading 
order of cumulants in traces of products of our compound Wishart 
matrices. For the usual Wishart matrices this was derived in [MN . 
Our more general version follows by the same kind of calculations (c.f. 
Capitaine and Casalis jCC, §5]). 

Theorem 3.5. Let {X^}j^ be a sequence of complex Gaussian random 
matrices. Put 

P N (Di) := X* N Df ] X N 
Let (k r ) r( zN denote the classical cumulants, then we have for all r G N 

(3) UmN^k^TrlPNiD^ ■ ■ ■ P N (D ni )], - ■ ■ , 

Tr[P N (D ni+ ... +nr _ 1+1 ) ■ ■ ■ Pjv(Aii+-+nJ]j 

= ^2 ifj^{d u ...,d ni+ ... +nr ) 

7reS W c(n(l),...,n(r)) 



Proof: Let n = n\ + • • • + n r and let 7 be the permutation with r cycles: 

(1, . . . ,ni)(ni + l, . . . ,n 1 + n 2 ) ■ ■ ■ («H hn r _i + l, . ..,nH \-n r ). 

Let 

Yi = Tr(P/v(Z) ni+ ..._|_„ i ■ ■ ■ P/v(D ril+ ... +ri; )) 
By jHTMl Theorem 2] 

(4) E(Yi ...Y r )=Y, N*^~ X ^~ n Tr a (Di, D 2 , . . . , D n ) 

For a 6 ^ let erV 7 be the partition of [n] whose blocks are the orbits 
of the group generated by a and 7. cr V 7 also defines a partition of the 
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r cycles of 7. Let us denote this partition of [r] by A 1 (a). Conversely 
let s be the number of cycles of a, a V 7 determines a partition of the 
cycles of a; we shall denote this by ^(7). Note that if A^(cr) = l r 
then a V 7 = l n and thus A a {pi) = l s . 

For a partition A = {Ai, . . . , A k } of [r] let 

k 

E x (n,...,y r ) = JjE(JJyi) 

If A = l r then E j4 (Yi, . . . , F r ) = E(Yi ■ ■ - Y r ). Equation (jlj) can now be 
extended easily to obtain 

E A (y 1 ,...,y r )= ]T iv^-^-^p!,...,^) 

A 7 (cr)< J 4 

Let S) be the Mobius function of the lattice of partitions; in 
particular /j,(A, l r ) = (-1) #(A) - 1 (#(A) - 1)!. Note that 

sail, 1° B < J ' 



= J] MAI^Ea^,...,^) 

AeV(r) 

= Yl MAlr) jV #(CT ~ l7) ~ nr Ir CT (.Di, ...,£>„) 

^ 7 (cr)<A 

o-eSn AeV{r) 

A-f{o)<A 

= N # ^ 1 ' y ^ n Tr (7 (Di, . . . , D n ) 

crV7=l„ 

= Y N*^- nr Tr t7 (D 1 ,...,D n ) 

crV7=l„ 

Recall that for er G 5 n with cr V 7 = l n there is an integer g = 
g{a) such that #(cx) + #(cr -1 7) + #(7) = n + 2(1 - g) and that a G 
S^c( n i, • • • 5 n r) means that cr V 7 = l n and g(a) = 0. 
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= N ~ n N * (a ' lj) N# (a hr r7 {D 1 ,...,D n ) 

creSn 
CTV7=1„ 

= N 2 - r -^hr a (D l} ...,D n ) 

<T£S n 

cV7=l„ 

= N 2 ^ J2 tr a (D u ...,D n ) + 0(N- r ) 

aeS N c(ni,...,n r ) 

Since linijvtr .(Di, . . . , D n ) = ip a (di, ■ ■ ■ , d n ) we have the required 
result. □ 

This theorem contains again all relevant information about the limit 
behaviour of the random variables T^P^-D^) • • • Pjv(-Dj„)). First, we 
have the following statement analogous to the law of large numbers. 

Corollary 3.6. The random variables {ir[Pjv(Z)j 1 ) • • • -P/v(-Di„)]}ii,...i„ 
converge in distribution to constant random variables ^(d^, . . . , di n ) ■ 1, 
where 

(5) (3(d h ,. . . ,d in ) := ^{d h ,...,d in ). 

Tr£NC(n) 

The form of ^(d^, . . . ,di n ) is, of course, in agreement with the fact 
that E{tr[Pjv(Ai) ' ' ' -P/v(A„)]} has, in the limit N — > oo, to agree with 
the corresponding moment of the compound free Poisson variables, 

ifj(c*d h c ■ c*d i2 c ■ ■ ■ c*d in c) 

where c is a circular random variable *-free from {d 1 , . . . , d p }. Again, 
we magnify the fluctuations around that limit, thus shifting the first 
cumulants to zero and getting only a non-vanishing limit for the second 
cumulants - hence getting normal limit fluctuations. 

Corollary 3.7. Consider the (magnified) fluctuations around the limit 
value, 

(6) 

F N (D h , . . . , AJ : = N- (tr[P N (D h ) ■ ■ ■ P N {D in )\ - [3(d h , . . .,dj) 

= Tr[P N {D n ) ■ ■ ■ P N (D in )\ - NP(d H ,. . . , d in ). 

The family of all fluctuations {F^{D il , . . . , D in )) ne ^ converges in dis- 
tribution towards a centered Gaussian family ( (F(di 1 , . . . , di n ) J , 
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with covariance given by 



(7) E{F(d h ,...,d im )-F(d, 




neS NC (m,n) 



Again, it remains to understand this covariance and we will be aiming 
at a more operator-algebraic description of these fluctuations in order 
to attack this combinatorial problem. 

4. Realization of semi-circular and free compound 
Poisson elements on Fock spaces 

The main theme of our investigations is the conviction that wher- 
ever planar or non-crossing objects arise, there is some free probability 
lurking behind the picture. Since the fluctuations of our Gaussian and 
Wishart random matrices can be described combinatorially in terms 
of non-crossing permutations, we expect also some operator-algebraic 
or some more abstract "free" description of this situation. Our main 
results in the coming sections will provide these descriptions. Let us 
begin by recalling the realization of a semi-circle and a compound Pois- 
son distribution on a full Fock space by using creation, annihilation, 
and preservation operators. 

4.1. Semi-circular case. For a real Hilbert space Hr with complex- 
ification H, we consider the full Fock space 



oc 



r(H) := 0H 8n = cn © n © n m © . . . 




and 



Kf) fl ® • • • © fn = f © fl © • • • © fn 

and the annihilation operator l*(f), by 

nm=o 



and 



r(/)/i©---©/„ = (/i,/)/ 2 ©---©/, 



n 



(neN,f 1 ,...,f n eH). 
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For / G Hr, we put 

u(f) :=/(/) + /*(/) 

and we will denote by A(Hr) the complex unital *-algebra generated 
by all u)(f) for / G Hr. Note that all u(f) are self-adjoint and that 
the vector fl is cyclic and separating for the algebra A(TCu). 
If we define on 7i an involution / i— > / by 



fi + ih ■= fi ~ if2 for f u f 2 G H R , 

then / i — > u(f) extends from a real linear mapping on TC R to a complex 
linear mapping on TC with 

uj(f) = i(f) + r(f) (fen). 

Note that the unital *-algebra generated by all ui(f) with / G TC is just 
•A(Wr). 

It is well known (see, e.g., VDNJ) that these operators u(f) have a 



semi-circular distribution and thus the asymptotics of the expectation 
of traces of Gaussian random matrices can also be stated as follows. 

Proposition 4.1. Let X^{f) (f G TCr) be a family of Hermitian 
Gaussian random matrices. Then for all fx, . . . , f n G T~Cr 

(8) lim E{tr[X N (fx)---X N (f n )]} = (u(fx)---uj(fn)n,ty. 

N^oo 

Let us in this context also recall the definition of the Wick products. 

Definition 4.2. For /i,...,/„6H the Wick product W(fx ® • • • <8> /„) 
is the unique element of A(TCu.) such 

(9) W(f\ ® • • • <g> / n )ft = /i <g> • • • ® /„. 

For n = 0, this has to be understood as VT(fi) = 1. 

Since Q is cyclic and separating for A(Hm), these Wick products 
exist and are uniquely determined. 

From the definition of the creation and annihilation operators it is 
clear that these Wick products satisfy for all f,fi,---,f n £ ^ the 
relation 

u(f)W(fx ® • • • ® f n ) = W(f ® A ® - • • ® f n ) + (fx, f)W(f 2 ® • • ■ ® f n ). 

This can also be used as a recursive definition for the Wick products 
and shows that W(fx ® • • • ® /„) is a polynomial in uj(fx), ■ ■ ■ , u(f n ). 

In the case f = fi = ■ ■ ■ = f n this reduces to the three-term recur- 
rence relation for the Chebyshev polynomials and shows that 

W(f*») = U n (u(f)/2), 

where U n is the n th Chebyshev polynomial of the second kind. 
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4.2. Compound Poisson case. In this case we start with a unital 
*-algebra T> equipped with a tracial state ip and represent T>, via the 

GNS-representation, on TC := T> ' \ where the inner product on TL is 
given by 

(di, d 2 ) := il)(d* 2 di). 
Then we take the full Fock space T{Tt) and consider there as before 
the creation and annihilation operators 1(d) and l*(d), respectively. But 
now we have, in addition, also to consider the preservation (or gauge) 
operator A(d) (d G T>) which is defined by 

A(d)Vl = 

and 

A(rf)/i ® • • • ® /„ = (c/fi) ® f 2 ® • • • ® f n 
for fi, . . . , f n G TC. Note that the multiplication DxP-tP extends 
to a module action 

VxH^H 
(d, f) ~ df. 

For deDwe define now 

(10) P {d) := /(d) + Z*(d*) + A(d) + ^(d)l, 

and we will denote by ^4(P) the unital *-algebra generated by all p(d) 
for d ET>. Note that we have 

= p(d*) for d E T>. 

One knows (see, e.g., |GSS| |SpH JNSp| ) that these operators p(d) 
give a realization of compound Poisson elements, i.e., their moments 
are given by Eq. ©. Thus we can state the asymptotics of the expected 
value of traces of our compound Wishart matrices also in the following 
form. 

Proposition 4.3. Suppose the family {d[ n \ . . . , Dp N ^} converges in 
distribution to {d ly . . . ,d p } in (V,ip). Let (X N ) N& be a sequence of 
N x N Hermitian Gaussian random matrices. Then 

(11) lim E{tr[P N (D h ) ■ ■ ■ P N (D in ))} = (jp{d h ) ■ ■ -p(d ln )n,n). 

Again, Wick products will play a role in this context. As before, 
these should be polynomials in the {p(d) \ d G T>} with the defining 
property that 

W(d 1 <g> ■ • • <8 d n )Vt = d x ® • • • <g) d n . 

However, in contrast to the semi-circular case, the multiplication of gPs 
in the arguments (under the action of A) has the effect that in order 
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to produce counter-terms for p(di) ■ ■ -p(d n )Q to get d± <S> ■ • • ® d n one 
also has to involve operators like p{d\d 2 ) etc. This means that W(d® n ) 
is in general not just a polynomial in p(d), but some polynomial in all 
{p(d k ) | k < n}. In particular, in general there is no relation between 
Wick polynomials W(d® n ) and the orthogonal polynomials with respect 
to the distribution of p(d). From the point of view of Levy processes 
the occurrence of p(d h ) in W(d® n ) is not very surprising, because this 
corresponds to the higher diagonal measures (Ito-formulas) and a Levy 
process should come along with its higher variations. 

It appears that Anshelevich pVnsJ was the first to introduce and in- 
vestigate these polynomials in this generality (and also some (/-deform- 
ations thereof). Since these polynomials appear implicitly in the clas- 
sical case in a paper of Kailath and Segall, he called them free Kailath- 
Segall polynomials. 

By taking into account the action of our operators on the full Fock 
space one sees quite easily that these Wick products should be defined 
as follows. 

Definition 4.4. For a given algebra T> with state ip, the Wick prod- 
ucts or free Kailath- Segall polynomials of the corresponding compound 
Poisson distribution are recursively defined by (d, d±, . . . , d n G T>) 

W(d) =p(d)-ip(d)l 

and 

W(d ® d l ® • • • ® d n ) = p(d)W(d 1 <g> d 2 ® • ■ ■ ® d n ) 

- ip(dd 1 )W(d 2 <g> d 3 (8) • • • ® d n ) 

- W(ddx g) d 2 <8> • • • <S> d n ) 

- ip(d)W(d 1 ®d 2 ®---®d n ) 

5. Cyclic Fock space 

Our main aim is to express the formulas for the limit fluctuations 
of Gaussian random matrices and of compound Wishart matrices also 
with the help of the operators uj(f) and p(d), respectively. In order to 
do so we have, however, to introduce another variant of a Fock space. 
Whereas the elements in the full Fock space, fx (8) • • • Cg) f n , are linear 
kind of objects - with a beginning and an end - we are looking on traces 
and thus should identify the beginning and the end in a cyclic way. 

Here are two versions of such a cyclic Fock space, the first one over 
arbitrary Hilbert spaces 7i and suited for semi-circular systems, and 
the second one over an algebra V and suited for compound Poisson 
systems. 
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Since for the calculation of our moments we only have to deal with 
elements in the algebraic Fock space (without taking a Hilbert space 
completion), we will restrict ourselves to this case in the following in 
order to avoid technicalities about unbounded operators. 

5.1. Semi-circular case. For a Hilbert space H, the algebraic full 
Fock space 

.Faig(w) := cn © n © n m © n m © • • • 

is generated by tensors /i©- • - ©/n, where we can think of the fi, ■ ■ ■ , f n 
as being arranged on a linear string. To stress this linear nature of the 
usual full Fock space, we will address it in the following as linear Fock 
space J-'imi'H)- In our tracial context, however, we should consider 
circular tensors, where we think of the fi, ■ ■ ■ , f n as being arranged 
around a circle. We will denote these circular tensors by [/i © • • • © f n ] 
and the corresponding n-th particle space by W®". If we pair two 
circles, then we have the freedom of rotating them against each other, 
so the canonical inner product for this situation is given as follows. 

Definition 5.1. The cyclic Fock space is the algebraic direct sum 

oo 

(12) **c(W) = 0?C 

n=0 

equipped with an inner product given by linear extension of 

(13) ([A © • • • © /„], [gi, © • • • © g m })cyc ■■= 

n-1 

8nm ■ ^(/l) 9l+k) ■ (f2, 92+k) ' ' ' (fn, 9n+k) , 



k=0 

where we count modulo n in the indices of g. 

Note that one can also embed the full Fock space into the cyclic one 
via 



1 

[fi © • • • © /„] = -f= Y] fk <E> fk+i ® • • • ® fi 



k-i- 



In order to write our formula for the fluctuations in terms of moments 
of operators we still use the operators on the linear Fock space, but we 
will in the end make things cyclic by mapping the linear Fock space 
onto the cyclic one. 

Definition 5.2. We consider the mapping c between linear and cyclic 
Fock space, 

c : jF lin (7i) — >■ JF cyc (7Y), 
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which is given recursively by 

(14) cfi = 0, c(f) = [f], 
and 

(15) c(/i<g>...<g>/ n ) = [/i®.. 



/n] + (/l,/n)-c(/ 2 



1!) 



fn- 



fa ® fa ® fa ® fa ® /s 



/l £ ft /J /* 

i i I I r 



Figure 2. 




[/i ® h ® h ® fa ® /s] (A, /s) [/ 2 ® /s ® / 4 ] (/i, / 5 ) (fa, fa) [fa] 

Figure 3 Figure 4 Figure 5 

Illustration of equation (15). Elements of the full Fock space, /i (8) 
/2®/3®/4®/5 are represented by linear "half pairings" (figure 2). The 
operator c takes a linear half pairing and wraps it around into a circle 
(figure 3). Then c pairs off the /'s until either one or none remains 
(figures 4 and 5) . The idea of a half pairing is a special case of a half 
permutation explained fully in [KMS]. 

Of course, one can also write down this explicitly, here are just two 
examples: 

(16) c (/i <g) h ® /a ® / 4 ® /s) = [/i ®fa®fa®fa® fa] 

+ (/i, / 5 ) • [/ 2 ® /s ® / 4 ] + / 5 ) • (fa, fa) ■ [fa] 

and 

(17) c (/i ® fa® fa® fi® fa® fa) = [fa ® fa® fa® fa® fa® fa] 

+ (fa, fa) ■ [fa ®fa®fa® fa] + (a, ; 6 > • (fa, fa) ■ [fa ® fa] 

Let us now consider the relation between this cyclic Fock space and 
fluctuations of Gaussian random matrices. So for the following, let 



20 



J. A. MINGO AND R. SPEICHER 



-^jv(J) be our Gaussian random matrices which converge, for N — > oo, 
in distribution to a semi-circular family, given by u(f) := 1(f) + l*{f) 
realized on the full Fock space. 

Our main point is now that we can express the fluctuations of the 
Gaussian matrices via the operators 00(f). 

Theorem 5.3. Let Xj^(f) (f e Hm.) be a family of Hermitian Gaussian 
random matrices. Then for all n,m e N and all fi, . . . , f n , gi, . . . , 
9m ^ T~Lr 

(18) lim hiTriX^h) ■ ■ -X N (f n )), Tr[X N (g x ) ■ -'X N {g m ))} 

iV^oo 

= (cu(fi) ■ ■■u(f n )Q,,cu(g m ) ■ ■ ■ uj(gi)ty cyc 

Note that the inversion of indices in the g's is forced upon us by the 
fact that our expression in random matrices is linear in both its traces, 
whereas our cyclic Fock space inner product is anti-linear in the second 
argument 

Remark 5.4. One might wonder whether the right-hand side of our Eq. 
(|18|) should not also have the structure of a variance. This is indeed 
the case, but is somehow hidden in our definition that cQ = 0. If, 
instead of c , we use the mapping c , given as follows 

ci] := ct] + (77, Q)Q, 

then the right-hand sided of (fTSj) has the form 

(C Tjt, C T] 2 ) eye = (c77l,C7? 2 } C yc ~ (c7/i,fi) cyc ■ (fi,C77 2 )cyc- 

We will prove Theorem 15.31 later as a corollary of a general limit 
theorem. For the moment, we will be content with checking the con- 
sistency of our statement with respect to traciality. Since the left hand 
side is tracial in the arguments of the traces, the right hand side should 
be tracial, too. Recall that A(Hk) is the unital *-algebra generated by 
all u(f) = l(f) + l*(f) with/ en. 

Lemma 5.5. The mapping 

(19) A(H R ) -> ^cyc(ft) 

(20) awcflfi 

is tracial, i.e., for all a,b e ^4(?^k) we have 

c abQ = c baVt. 

Proof. Since fl is cyclic and separating for the algebra A(TCm), it suffices 
to show that 

cu(f)W(f 1 <g> • • • (8) f n )Q = c W(h ® ■ ■ ■ ® f n )uj(f)Q 
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for all f,fi,---,f n £ On the left side, we have 

CW(/)^(A® - ■ -®/„)n = C (/®/!® - • -®/n) + (/l, /)'C(/ 2 ®. •■»/„). 

For the right side, it follows from the relation 1 <E> ■ ■ ■ <8> / n )* = 

W(f n ® • • • <g> /i) that 

^®---®/«M/) = m®---®/ n ®/) + (/„JW/i8---®U, 

which yields 

C • -®f n Mm = C (/l®- • -®f n ®f) + (fn, /)-C (A®" • -®/n-l)- 

From the definition of c we see that both sides are the same. □ 

5.2. Compound Poisson case. Let us now consider the case where 
we have a *-algebra V with trace ip. We denote by 

oo 

^iinOP) = 02^ n 

n=0 

the algebraic linear Fock space and by 

oo 
n=0 

the algebraic cyclic Fock space. 

Since in this case we also have actions of our operators which multiply 
inside the argument, we have to take this into account when we glue 
the beginning and end of the tensors together. Thus we have to change 
the definition of the map c as follows. 

Definition 5.6. We consider the linear mapping 

c : Fn n (V) -> F cyc {V), 

given recursively by 

cfi := and c(d) = [d] 

As we have been unable to locate a proof of this in the literature one is given in 
KMS §10]. The main idea is to write w(/i)w(/2) • • -u>(f n ) as a linear combination 
of Wick polynomials. In fact we may write w(/i)cj(/ 2 ) ■ ■ ■ w(/„) = J2ir W n (f\® fz® 
■ ■ ■<£>/„) (*) where the sum is over all non-crossing "half-pairings" of [n], that is non- 
crossing pairings of [n] with only singletons and pairs such that the singletons are 
not enclosed by any pair, and W w (/i®-- -®/ n ) = (/in/ji) ■ • • (fi k , fj k )W(f h ®- ■ -(g) 
fl m ) where the pairs of ir are (ii, ji ),..., (ik, jk) and the singletons are . . . , (l m ). 
Note that in (*) there is only one term involving a Wick polynomial on n fs. So 
that if we know that W(fi ® • • • <8> f m )* — W(f m ® • • • ® /i) for m < n they we 
apply induction to the adjoint of (*). 
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and 

(21) c (di ® • • • <g> d n ) :— [dx <g) ■ • ■ <g> d n ] + [d n d\ ®>d 2 ® ■ - -® d„_i] 

+ ip(did n ) ■ c (d 2 <S> ■ ■ ■ <8> 

For example 

(22) c (g?! ®d 2 ® d 3 ) = [di ®d 2 ® d 3 ] + [d 3 d 2 ® d 2 ] + [<fc] 




[di <®d 2 ® d 3 ) [d 3 di ® d 2 ] ip{d 3 di) [d 2 ] 



Figure 7. Figure 8. Figure 9. 

Illustration of Equation (|22|) The vector di®d 2 ® d 3 is represented 
by a linear half permutation (figure 6) with one open block for each 
factor in the tensor product (c.f. lKMS.1 §10]). [d± <8> d 2 <g> d 3 ] is rep- 
resented by a circular half permutation with one open block for each 
factor in the tensor product. ip(d 3 d\)[d 2 ] is represented by a circular 
half permutation with one closed block and one open block (figure 9). 
The operator c first turns the linear half permutation into a circular 
half permutation (figure 7). Then c fuses a pair of open blocks (figure 
8) and then closes the just formed open block (figure 9). This process 
continues until either one or zero open blocks remain. 



Then we claim that one can express the fluctuations of our compound 
Wishart matrices also by calculations in terms of the corresponding 
operators p(d). 

Theorem 5.7. Suppose {D[ N \ . . . } Dp N ^} converges in distribution to 
{d\, . . . , d p } in (V, ip) and (X^Nm is a sequence of N x N Hermitian 
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Gaussian random matrices. We put P/v(A) := X^D^ and let 

p(di) be our operators on the full Fock space, then we have for all m,n G 
N that 

(23) lim h{Tr[P N {D h ) ■ ■ ■ P N (D ln )], Tr[P N (D ln+1 ) ■ ■ ■ P N (D ln+ J]} 

iV— +00 

= (cp(d h ) ■ ■ -p(d in )n, cp(d in+1 ) ■ ■ •p(rfi n+m )^)cyc- 

Again, we check only the traciality of the right hand side and post- 
pone the proof of the statement until we have proved our general limit 
theorem. Recall that we denote by A(T>) the unital *-algebra generated 
by all p(d) for d G £>. 

Lemma 5.8. The mapping 

A(V) -> F cyc (V) 
a 1 — ^ c aVL 

is tracial. 

Proof. Since fl is cyclic and separating for A(V) (see jAnsj ). it suffices 
to check for d, d\, . . . , d n G T> that 

cp(d)W(di <g> • • • ® d n )VL = c W(di <8> ■ ■ ■ <S) d n )p(d)Vt. 

For n = 0, i.e., W(Q) = 1, this is surely true. In general we have for 
the left hand side 

cp(d)W(d 1 <8 • • • <8 d n )VL = cp(d)dt <8> • • • <8> d n 

= c (d <S> di ® ■ ■ • <g> d n + ip(ddi)d 2 <8> • • • (8 d n 

+ ddi ® d 2 <8> • ■ ■ <8> d n + ip(d)d 1 (g) • - • ® d n ) 

By using the identity 2 W(di (8> • • • ® d n )* = VF(d* (8> • • • <8 d*) we have 

W(di <8 • • • <8> d„)p(d) = (8 • • • <8> d n <8 d) 

+ iy (d x ® • • • (8) d n _i (8) d n d) 
+ V>(d)W(di® 

Thus the right hand side becomes 

c W(d 1 <g> ■ ■ ■ (8) d n )p(d)Vt = c (d x <8 • - • <8> d n ® d 

+ ijj(d n d)d 1 (8) • • - <8> dn-i 
+ di (8) • • • <8> d„_i (8 d n d 
+ ^(d)di ® • • • <8 d n ) 

2 The proof is very similar to that sketched in the footnote on page|2J A detailed 
proof is provided in §10]. 
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So it remains to show that 

c (d ® di <g> • • • ® d n + il)(ddi)d 2 <g> • • • ® d n + ddi ® d 2 ® ■ ■ ■ ® d n ) 
— c (di ® • • • ® d n ® d + ip(d n d)d 1 ® • • • ® d n _ x + d x ® ■ ■ ■ ® d n _i ® d n d) 

This can be checked directly by applying the definition of the mapping 
c. □ 

6. Second order freeness and abstract limit theorems 

We shall derive our main theorems. 15. HI and 15. 7\ from a general limit 
theorem, very much in the same spirit as one can get the distribution 
of the semi-circle and the compound free Poisson distributions from 
free limit theorems, see |Spl| . The crucial idea is the notion of second 
order freeness which we introduce below. 

Definition 6.1. A second order non- commutative probability space (A, 
tp, p) consists of a unital algebra A, a tracial linear functional 

ip : A -> C with ip(l) = 1 

and a bilinear functional 

p : A x A -> C, 
which is tracial in both arguments and which satisfies 
p(a, 1) = = p(l, b) for all a,b G A. 

Notation 6.2. Let unital subalgebras Ai, . . . , A r C A be given. 

1) We say that a tuple (ai,...,a n ) {n > 1) of elements from A is 
cyclically alternating if, for each k, we have an i{k) e {1, . . . , r} such 
that at G Ai{k) and, if n > 2, we have i(k) ^ i(k+l) for all k = 1, . . . , n. 
We count indices in a cyclic way modulo n, i.e., for k = n the above 
means i(n) ^ z(l). Note that for n = 1 we mean that a\ is in some Ai- 

2) We say that a tuple (ai, . . . , a n ) of elements from A is centered if we 
have 

¥{ a k) — for all k = 1, . . . , n. 

Definition 6.3. Let (^4, <p, p) be a second order non-commutative prob- 
ability space. We say that unital subalgebras Ax, . . . , A r C A are free 
with respect to {ip, p) or free of second order, if they are free with re- 
spect to ip and whenever we have centered and cyclically alternating 
tuples (di, . . . , a n ) and . . . , b m ) from A then we have: 

i) p(ai ■ ■ ■ a n , b\ • • • b m ) = for n ^ m; 
ii) p(a, b) = for a G Ai, b G Aj, and i ^ j; 
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Hi) if n = m > 2, then 



n-l 



p(ai • • • a n , bi 



fc=0 



Note that in the sum the indices of the a's increase, whereas those 
of the 6's decrease; one should think of two concentric circles with the 
a's on one of them and the 6's on the other. However, whereas on one 
circle we have a clockwise orientation of the points, on the other circle 
the orientation is counter-clockwise. Thus, in order to match up these 
points modulo a rotation of the circles, we have to pair the indices as 
in the sum above. 

Condition (Hi) is the annular version of the disc picture of first order 
freeness: suppose a& G Ai k with 4>(a,k) = for 1 < k < n and we ar- 
range the elements around the boundary of a disc (figure 10). 
The only non-crossing partition of [n] that only connects elements from 
the same algebra consists of all singletons and since of a singleton is 
0, we have that <f>(ai ■ ■ ■ a n ) = 0. 

In the annular case we put the centered and cyclically alternating 
elements ( ) and (pi, ... , b n ) around the boundary of an an- 

nulus. We only connect elements from the same algebra and as the 
elements are centered we have no singletons; so we must connect in 
pairs elements from opposite circles is all possible ways (figure 11). 
This is the meaning of condition (Hi). 



Note that, as in the case of freeness, the trick of writing elements a 

as 



allows us to calculate p in terms of cp and p restricted to the subalgebras. 
However, whereas the formulas for tp of mixed moments contain only 
ip applied to the subalgebras, p of mixed moments has in general to be 
expressed in both ip and p restricted to the subalgebras. 




Figure 10. 



Figure 11. 



a = a° + ip(a) ■ 1, 



where y?(a°) = 0, 
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For example, assume we have two subalgebras Ai and A.%, and ele- 
ments ai,a 2 G A\ and b±, b 2 G A 2 . Then we have 

(24) p (a l ,b 1 ) = 0, 



(25) p(a 1 b 2 ,a 2 ) = ip(b 2 ) ■ p(a 1 ,a 2 ) 
or 

(26) p(ai&2, a 2 b 2 ) = ip(a l a 2 )ip{b 1 b 2 ) - ip(a 1 a 2 )ip{bi)ip{b 2 ) 

- ip(ai)(f(a 2 )ip(b 1 b 2 ) + ip(a 1 )ip(a 2 )(p(bi)(p(b 2 ) 

+ p{a u a 2 )(p(bi)cp(b 2 ) + cp{ai)(p(a 2 )p{b u b 2 ). 

One should note that each variable appearing in the arguments of 
p on the left-hand side of these examples has to appear exactly once 
in each product on the right-hand side. Let us formalize this in the 
following definition. 

Notation 6.4. Let (A, (p, p) be a second order non- commutative prob- 
ability space with subalgebras Ai, . . . , A r C A, and consider elements 
ai,...,a n G L)l =1 Ai- A balanced expression (with respect to the sub- 
algebras A\, . . . , A r ) in ax, . . . ,a n is a product of factors <p(ai 1 ■ ■ ■ a it ) 
and p(oj 1 • • • aj s , aj 1 ■ ■ ■ aj t ) where each aj has to appear exactly once 
among all arguments and the argument of each tp or the arguments of 
each p contains only Oj's from a single Aj. 

For example, balanced expressions in ai, a 2 , as, are 
ip(aia^)(p(a 2 a4) if 01,03 G A± and a 2 , a 4 G A2 
or ip(ai)tp(a4)p(a 2 , a 3 ) if ai, a 4 G Ai and a 2 , 03 G A2 
Every summand on the right-hand side of Eq. (|2Tj) is a balanced ex- 
pression in ai, a 2 , 61, 62 if ai, a 2 G Ai and 61, 62 G A- 

Lemma 6.5. Let Ai, . . . , A r in (A, ip, p) be free with respect to (</?, p). 
Suppose we have cyclically alternating (ai, . . . , a n ) and (61, ... , b m ) and 
denote by s the number of different subalgebras appearing in {a\, . . . a n , 
bi, . . . b m }. Then p(a± ■ ■ ■ a n , b\ ■ ■ ■ b m ) is either or can be written as a 
sum of balanced expressions in a\, . . . , a n , 61, . . . , b m , such that each of 
these balanced expressions has at least s factors and contains at most 
one p-factor. 

Thus any expression of the form p(ai ■ ■ ■ a n , bi ■ ■ ■ b m ) for cyclically 
alternating (ai, . . . , a n ) and (61, ... , b m ) is determined by the value of 
p restricted to A\ U ■ ■ • U A r and by the value of p restricted to (A x 

A) u • • • u (A x AO- 



SECOND ORDER FREENESS 



27 



Proof. We will prove this by induction on n + m. The case n + m = 2, 
i.e., n = m — 1, is clear. 

So consider n + m > 3. Put 

a£ := a fc - <p(a fc ) ■ 1, b° := b x - <p(bi) ■ 1 

for k = 1, . . . , n and I = 1, . . . ,m. Then we have 

p(ai ■■■a n ,b 1 -- - b m ) 

= p(K + ^( ai )) • • • « + <^(a n )), (6? + ^(60) • • • (b° m + ip(b rn ))) 
= ^K(i)) • ■ ■ <p(<hw) ■ f(b q(1) ) 

(27) 

X • • ■ X <P{bq(l)) ■ P{a°p(l) ■ ■ ■ Op(n-fc): • • ■ & 5(m-i))' 

where the sum runs over all partitions 

((p(l), • • -iP( k ))i (Pi 1 ), ...,p(n- k)) of the set [l,n] 

and 

((g(l), . . . , g(0), (g(l), ■ ■ ■ , q(m - /)) of the set [1, m\ 

into two ordered subsets (with k = 0, . . . , n and I = 0, . . . , m). The 
term corresponding to k = I = is, by Definition IO| either (when 
m ytz n) or is a balanced expression in the centered elements with 
at least one factor for each occurring subalgebra. Now note that a 
balanced expression in centered elements can be rewritten as a sum of 
balanced expressions in the original elements and that the number of 
factors can only increase by doing so. 

For the other terms with k + l > 1, • • • , a p( n -k)) anci (^(i)> • • • > 

^/ TO _n) may no longer be cyclically alternating. So we group together 
adjacent elements from the same algebra to produce a cyclically alter- 
nating tuple with at least max{l, s — (k + I)} subalgebras appearing, 
and so we can apply our induction hypothesis. Indeed, the term 

(28) P( a p(l) ' ' ' a p(n-k)y b°q(i) ■ ■ ■ ^( m ^)) 

contains elements from at least s — (k + l) different subalgebras; by our 
induction hypothesis, we may write it as the sum of balanced expres- 
sions in the a°'s and b 0, s, each product containing at least s — (k + I) 
factors. Again we write a balanced expression in centered elements as 
a sum of balanced expressions in the original elements. This means we 
can write the term (}2"%|) as a sum of balanced expressions in 

a p(l)i " " " j Qp(n-fc), ' ' ' j bq(m-l) 
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with at least s — (k + I) factors for each product. Together with the 
k + / factors 

^K(i)) • • • ¥>(op(*)) • v(6 9 (i)) • • • ¥>(&«(!)) 

this gives the assertion. Note that in all our steps balancedness is 
preserved and that at most one p-term can occur in all the reductions. 

□ 

A very special case of such a factorization is given in the next lemma. 

Lemma 6.6. Let (A, <p, p) be a second order non- commutative probabil- 
ity space and let A±, . . . , A r C A be free with respect to (if, p) . Consider 
cyclically alternating (oi . . . , a n ) and (bi, . . . , b m ) from A. 

1 ) Assume that the subalgebra of ai appears only once. Then we have 

p(ai ■ ■ ■ a n , h ■ ■ ■ b m ) = ip(a 1 )p(a 2 ■ ■ ■ a n , h ■ ■ ■ b m ). 

2) Assume that the subalgebra of b\ appears only once. Then we have 

p(at ■ ■ ■ a n , h ■ ■ ■ b m ) = <p(bi)p(ai ■■■a n ,b 2 --- b m ). 
Proof. We only prove the first part. Put 

a° k := a k - <p(a k )l, 6° := k - <p(b t )l. 

We have 

p(a 1 ■ ■ -a,,, h ■ • • b m ) = p((a° + y?(ai)l)a 2 • • • a n , 6 X • • • b m ) 

= p(a°a 2 • • • On, h • ■ ■ b m ) + ip(a 1 )p(a 2 ■ ■ ■ a n , b x ■ ■ ■ b m ). 

We shall show that the first term is 0. 

Indeed, we shall show that if (ai, a 2 , . . . , a n ) and(6i, b 2 , ■ ■ ■ , b m ) are 
cyclically alternating and the algebra of a\ appears only once then 
p(a°a 2 ■ ■ • a n ,bi ■ ■ • b m ) = 0. We shall do this by induction on m + n. 
By equations and we have p(a°a 2 ■ ■ ■ a n , b± • • • b m ) = when 
m + n = 2 or 3. Suppose we have proved the result for m + n < j; we 
shall prove it for m + n = j . 

We shall use the expansion in equation (|27|) and show that for 1 < 
k < n - 1 and 1 < I < m, p{a°a°^ ■ ■ ■ a° (fc) , 6° (1) • • • b°^) = for all 

subsets {p(l), • • • ,p(k)} C {1, 2, 3, . . . , n - 1} and {q(l), . . -,q(l)} C 
{1,2,3,. ..,m}. 

When k = n — 1 and / = m we have that (a°, . . . , a°) and (6°, . . . , b^A 
are centered and cyclically alternating. If m ^ n we have p(a°a2 ■ ■ • 
a°, b° ■ ■ ■ b^A = by (i) of Definition 16.31 If m — n then by (Hi), we 
have p(a°ct2 ' ' ' fl ni b°-- - 6^) = because <p(a°6°) = for all i. 

Suppose next that k + I < m + n — 2. We can no longer expect 
K> a p(i)> ■ ■ ■ > a p(fc)) and • • • > & g(o) to be cyclically alternating; so 



SECOND ORDER FREENESS 



20 



we group adjacent terms from the same algebra and write a°a°^ ■ ■ ■ 
a ° P {k) = a ° c i " ' °s and b q(i) ■ ■ ■ b ° q (i) = d x ---d t with (ci, . . . , c s ) and 
(di, . . . , d t ) cyclically alternating and neither involving the algebra of 
a\. Now s + t + l < k + l + 1 < m + n — 1 < j — 1, so by our induction 
hypothesis p«a° (1) ■ • • a° (fc) , b° q{1) ■ ■ ■ b° q{l) ) = p(a°ci • • - c s , d x ■ ■ ■ d t ) = 

□ 

If by the lemma above we successively remove all subalgebras which 
occur only once and multiply together cyclic neighbours from the same 
subalgebra, then we arrive finally at p(a, b) for a, b from one of the 
subalgebras (both from the same, in order to get a non-vanishing con- 
tribution) or at p(a± • • • a n , b± • • • b rn ) where both arguments are cycli- 
cally alternating and in addition each involved subalgebra appears at 
least twice. In the latter case we have either a very special matching 
of the involved subalgebras or we can strengthen Lemma fo. 51 to obtain 
at least one more ip-f actor. 

Lemma 6.7. Let (A, <p, p) be a second order non- commutative probabil- 
ity space and let A\ , . . . , A r C A be free with respect to (ip, p) . Suppose 
(ax, ... , a n ) and (&i, . . . , b m ) are cyclically alternating and denote by s 
the number of different subalgebras appearing in {ai, . . . , a n , b%, . . . b m }. 
Suppose also that each involved subalgebra appears at least twice. 

Then p(a\ ■ ■ ■ a n , b\ ■ ■ ■ b m ) can be written as a sum of balanced ex- 
pressions with at least s + 1 factors unless the following conditions are 
satisfied: 

— m = n = s; 

- for each k there is k' such that au and bk> are from the same 
subalgebra; and 

there is q such that for all k, k' = —k + q mod n. 

In this case p{a\ ■ • • a n , b\ • ■ ■ b n ) = (p(a\by) ■ ■ ■ ip(a n b n >) + S, where S is 
a sum of balanced expressions with at least s + 1 factors. 

Proof. Let us look again at the expansion 

p(ai ■■■a n ,bi-- - bm) 

p,q 

(29) ■ p(ap(i) ■ ■ ■ dp(n-k)i ^(i) ' ' ' bq(m-i)) ■ 

First, consider a term with k+l > 1. Then there are two possibilities. 
If all {a p (i), . . . , Op(fc), • • • , b q (i)} belong to different subalgebras, 
then there must be exactly s subalgebras in {eip(i), . . . , %( n _fc), &?(i), 
} because each involved subalgebra appears at least twice. 
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If we group together any adjacent terms that may come from the same 
subalgebra we obtain cyclically alternating arguments and so can apply 
Lemma IHISl According to Lemma Id^I we can write p(a^ 1 - ) • • •fl& n _fc\, 
6°^ • • • &|( m _n) as a sum of balanced expressions with at least s factors. 
Combining these with the k + I factors <p(a p m) ■ ■ ■ <f(a p ^))<f(bq(i^) ■ ■ ■ 
(p(b q n\) we have that every term with k + l > 1 can be written as a sum 
of balanced expressions with at least s + 1 factors. 

Second, consider the term p(a° • • • a° , 6° • - • fo^J corresponding to k — 
1 = 0. If m 7^ n we have this is zero by Definition l6.31 So suppose m = 
n. Again by Definition O P« •■<,&?"• K) = EVo fKK+k) ' ' " 

If s < n then each term <p( a i^n+fc) ' ' ' f( a nH+i) nas n> s + 1 factors 
or is zero if for some factor <p(a°6° +fc _ r ), a r and 6 n +fe-r come from 
different algebras. Thus we get either or a balanced expression with 
at least s + 1 factors. 

Finally assume that s = m = n. Each subalgebra must appear 
exactly twice, so for each k there is k' such that at and by are from the 
same subalgebra, or else for all q, tp(a\bn +q _i) • ■ -ip(a°b°) = 0. Again 
we will have <p(°i^n+g-i) ' ' ' v( a nbq) — unless k' = —k + q mod n, for 
some q. For this q we have p(a° ■ • • a°, b° ■ • • 6°) = <p(a°b°,) ■ ■ ■ <p(a°6°,). 

By substituting <p(a^6^) = ip(akbk>) — <p(a,k)<p(bk') into ip(a°b°,) • • • 
P( a nK') we ma y write P( a ° • • • Or) &i • • • b°) as <p(aifoi') • • • <p(a n 6 n /) plus 
a sum of balanced expressions with at least n + 1 = s + 1 factors. □ 

We are now almost ready for the main limit theorem of second order 
freeness. It will turn out that moments of the limit can be calculated in 
terms of annular non-crossing objects. However, in this setting we will 
not arrive directly at permutations (as in the fluctuation formulas for 
random matrices), but - as is much more natural in the context of limit 
theorems - at partitions. In the random matrix setting of section 3 we 
got contributions of the form ip n for non-crossing permutations 7T. So 
we have to define the analogous object ip a for non-crossing partitions 
a. However, for non-crossing partitions, the contribution to ip a of a 
block which is the only through-block will require special treatment. 
We will need two different types of functions, ipi in the case of multiple 
through-blocks and tp 2 in the case of a single through-block. 

Notation 6.8. Let T be an index set and let two functions 
fa : |Jt"^c, (h, . . . ,t n ) ^ fa(ti, . . . ,t n ) 

neN 
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and 

^ 2 : [J T n xT rn -*C 

n,mefi 

(ti, . . . , t n ) X (t n+ i, . . . , t n+m ) I— > 1p2\tlj ■ ■ ■ -i tn'i tn+ly • • • ; tn+m) 

be given. Assume that ipi is tracial in its arguments, i.e., for all n G N 
and all t 1: . . . , t n G T we have 

i>l(tl,t 2 , ...,£„)= ^l(t 2 , • • • ,tn,h), 

and that ^ 2 is tracial in each of its groups of arguments, i.e., for all 
n, m G N and allti, . . . , t n+m we have 

ip2\ti, ■ ■ ■ , t n ; t n+ i, . . . , t n+m ) = ip 2 (t 2 , ■ ■ ■ ,t n , h] t n+ i, . . . 

and 

^2(^1, ■ ■ ■ ,t n ; t n+ i, . . . , t ra+rrt ) = ip2{ti, ■ ■ ■ ,t n ; t n+ 2, • • • , t n+m , 

Fix n, m G N and consider an annular non-crossing partition <x G 
NC(n,m). Then, for given t±, . . . , t n , t n+ ±, . . . t n+m G T we define 
ip a (ti, . . . ,t n ] t n+ i, . . . , t n+m ) as follows: If is not the only through- 
block of a then we choose the unique cyclic order on B (c.f. section 
2.2) and, writing it as a cycle B = (i(l), . . . ,i(k)), we put 

(30) ipB(ti, ■ ■ ■ , t n ; t n+ i, . . . , t n+m ) := tpi (ti(i),ti(2), ■ ■ ■ , ti(k))- 

If B is the only through-block of a, then we write it as B = BiVJB 2 with 
B 1 = (i(l),...,i(fc)) C [l,n] and B 2 = . . . C [n+l,n + m], 

where we induce the cyclic order of [l,n] on B\ and the cyclic order of 
[n + 1 , n + m] on B 2 . For such a block 5 we put 

(31) ipB(ti, ■ ■ ■ ,t n ;t n+ i, . . . ,t n+m ) := ^2(^(1), • • • ,ti(k)',tj(i), ■ ■ ■ ,tj(i)). 

Finally, we define 
(32) 

V'it(^I) • • • 5 t n , t n _|_i, . • . , tn+m,) . J^J '0_b(^1) • • • ; ^nj ^ra+l; • • • ; tn+rri)- 

Here are some examples of this notation: Consider n = m = 2 and 

a, = {(1,3), (2, 4)} (7 2 = {(1,2, 3), (4)}. 
Then <7i has two through-blocks so 

(33) ^ CT1 (ti,t 2 ;t 3 ,t 4 ) = t 3 )ipi{t 2 , U) 
whereas <r 2 has one through-block so 

(34) 4>a 2 (ti,t 2 ;t 3 ,t4) = ^i(t4)^ 2 (h,t 2 ;t 3 ). 
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Theorem 6.9. Let (An,<Pn, Pn) (N G N) be second order non-com- 
mutative probability spaces and let, for each N G N, unital subalgebras 
A N , . . . , A% C An be given which are free with respect to (<Pn, Pn)- Let 
T be an index set and assume that we have, for each t G T and each 
iVeN, elements 

ct N (t)eA N (i = l,...,JV), 
such that the following properties are satisfied: 

(a) The distribution of the q N (t) under ((Pn,Pn) is invariant under 
permutations of the upper indices, i.e., for all N G N, and all permu- 
tations n : [1, N] — > [1, N] we have for all n, m G N, ii, . . . , t n+m G T 
and all ■ ■ ■ , i(n + m) G [1, N] that 

Vjv(?a t (h) ■ ■■q l N n \t n )) = <^jv(g^° i(1) (ti) ■ ■■qT {n \ t n)) 

and 

PN(q N l \ti) ■ ■ ■ q N n \t n ) , q N n+l \t n+ i) ■ ■ ■ q N n+m \t n+m )) 

PN [Qn { t l)"'QN {tn),qN v ) ' ' ' Qn {tn+m)) 

(b) For all n, m G N and all t\, . . . , t n , t n +\, . . . , t n+m G T there exist 
constants ipi(t 1 , . . . ,t n ) and tp2{ti, ■ ■ ■ ,t n ;t n+ i, . . . ,t n+m ) such that 

(35) lim N ■ i p N (q N (t 1 )---q i N (t n )) =Mt u ...,t n ) 

N — >oo 

and 

(36) lim N ■ pjvfe(*i) " " " ?M*n+i) • ■•?»(*»+»)) 

TV— >oo 

— V^^lj ■ ■ ■ , ^nj ^n+lj ■ ■ ■ j tn+m). 

Fort eT and N EN let 

S N (t) := g^(t) + ■ • • + g#(t) G A N . 

Then we have 

(37) lim p N (S N (ti) ■ ■ -S N (t n ), S N (t n+1 ) ■ ■ ■ S N (t n+m )) 

iv— >oo 

^ ■ ■ ■ i^ni tn+li • • • j tn+m) 

a<=NC(n,m) 

Note that the left-hand side of the expressions (|33j) and (JSfiJl are 
independent of the value of the index i, and that the functions tpi and 
ip2 defined there have the traciality properties which are required in 
Notation Ol 
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Proof. For better legibility, we will suppress in the following the index 
iV at ifN and pjq and just write cp and p, respectively. 
We have 

E p{qf\h)---qi n \t n ),q^ n+1 \t n+1 )---q t £ ! n+rn \t n+m )) 
i:[l,n+m]-+[l,N] 

Because of our invariance assumption (a), the value of the term 

(38) p(qf\ tl ) ■ ■ ■ qf\t n ), qf +1 \t n+1 ) ■ ■ ■ qf +m \t n+m) ) 

depends on i only through the information where these indices are 
the same and where they are different. As usual, this information is 
encoded in a partition a of the set [1, n+m], and we denote the common 
value of for all i with ker(z) = a by 

(39) Pa{qN(h) ■ ■ ■qN(t n ),QN(t n +i) ■ • • qN(t n +m)) ■ 
Then we can continue our calculation as follows: 

= E E p{4 l \h)---qf\tn\ 

a&V(n+m) i:[l,n+m]— >[1,N] 

ker(i)— a i(n+l) /, \ i(n+m) i , \\ 

IN \ l n+l) ■ ■ ■ qN { l n+m)) 

= E P<r(lN{tl) ■ ••?Ar(*n),?Ar(*n+l) " " -9iv(Wm)) • {N)\ a \, 
<r£'P(n+m) 

because the number of i : [1, n + m] — > [1,N] with the property ker(i) = 
cr is given by 

N(N-l)---(N-\a\ + l)=: (JV) H . 

We have now to examine the contributions for different a. Let us 
first assume that cr contains a block B which is either contained in 
[1, n] or contained in [n + 1, . . . , n + m] and all of whose elements are 
consecutive in the induced cyclic order. Because of traciality of p it 
suffices to consider the case B = [1, s] for some s with 1 < s < n. By 
Lemma f6. 61 this implies 

P<r(qN(tl) ■ ■ ■ <lN(t s ) ■ ■ ■ C/iv(t n ), C/jv(t n+1 ) ■ • • C/Ar(t n+m )) 

= if{q N {t 1 ) ■ ■ ■ q N {t s )) ■ p a , (q N (t s+1 ) ■ ■ ■ q N {t n ), q N {t n+1 ) ■ ■ ■ q N (t n+m )) , 

where a' is that partition which results from a by removing the block 
B = [1, s] and relabelling elements. Since 

lim N ■ (p(q N (ti) ■ ■ ■ q N (t s )) = • • -,t s ), 

TV— >oo 
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the block B makes exactly the contribution to the final result as claimed 
in Eq. (|37|). Thus, by successively removing such blocks, it suffices to 
consider a's which have no blocks which are contained in either [l,n] 
or [n + 1, n + m] and which consist of cyclically consecutive elements. 

So let us now assume that a contains no blocks which are contained in 
either [1, n] or [n + l,n + m] and which consist of cyclically consecutive 
elements, and consider (|39|). By multiplying together neighbouring 
elements corresponding to the same block of a we can rewrite the two 
arguments of p in a cyclically alternating form. The fact that a contains 
no blocks of the form treated above implies that after this rewriting 
of arguments each involved subalgebra occurs at least twice. But then 
Lemma f6.7l implies that, unless condition (*) is satisfied, we can write 
all these terms as sums of products of at least |<r| + 1 factors. By our 
assumption, each of these factors multiplied by N converges to a finite 
number; however, since we have more than \a\ factors, this product 
multiplied by iV'""' will vanish in the limit N — > oo. This means that 
we can only get a non-vanishing limit for a a which satisfies condition 
(*) of Lemma 16.71 However, these are exactly the cases where each 
block B of a is of the form B = B\ U B2, where B\ C [l,n] and 
B2 C [n + 1, m + 1], are non-empty, and each consists of consecutive 
numbers with respect to the inherited order. Furthermore, the cyclic 
order of the restrictions of all blocks to the interval [1, n] must be 
the inversion of the cyclic order of the restrictions of all blocks to the 
interval [n + 1, n + m]. In this case (I39|) calculates as follows. If we 
have only one block in a, then our assumption, Equation ((HEl), gives, 
in the limit, for such a a the contribution 

VM^Ij ■ ■ ■ , t n ; t n+ i, . . . , t n+m ). 

If, on the other side, a has more than one block, then we get, according 
to the description of annular non-crossing partitions in section 2.2 and 
our assumption the product of ipB over all blocks B of a. 

Note that the reduction above leads to non-vanishing contributions 
exactly for non-crossing partitions a from NC (n, m) and each such 
partition a contributes a term ip a (ti, . . . ,t n ; t n+ i, . . . , t n+m ). □ 

7. Proofs of Theorems 15.31 and 15.71 

Now we can prove our main theorems by reducing them to the situ- 
ation covered in our limit theorem. 
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7.1. Proof of Theorem 15. 3L We have to show that for all n, m G N 

and fx , . . . , fn+m £ 

(cu(fi) ■ ■ -uj(f n )Q,cuj(f n+m ) ■ ■ -U (f n+ i)ft) cy C 

= e n (/*./i>- 

■K&NC2 (n,m) (i,i)E7r 

Note that we can, for any iVeN, replace 7i by 

A' 

ft = ft © ... © ft (JV summands) 
i=i 

and by 

1 w(/e ••■©/)• 



We can then put this into the framework of our general limit theorem 
by letting T = ft K , 

N 

A N = A{Q)H), A* n = A(Q@---@H@---@0) 

i=l 

<Pn{cl) = (aCl,Q) (a G .Ajv) 

and 

PAr(a, 6) = (coil, cfr*f2) cyc (a, 6G.Ajv) 
and finally, for / G ft K , 

<&(/) = J=u;(0 9-©/®-0O)g4 

ViV i-th 



Let us check that A N , . . . , A% C An are free with respect to (<Pn, Pn)' 
Freeness with respect to tp n is well-known, so we only have to consider 
p N . Take centered and cyclically alternating tuples ( Oi, • • • , CL n ) and 
(6i, . . . , b m ) from .4.^- Let us only consider the case n, m > 2, the cases 
were at least one of them is 1 are similar. Note that the centeredness of 
the aj implies that each a^l has no component in the direction Q and 
thus, by the fact that neighbours are from algebras with orthogonal 
Hilbert spaces, we have 

aia2 • ■ ■ a n Q = (aifl) © (a^fi) Cg) • • • <8> (a n Q). 

Since also the first and the last element are orthogonal, the action of 
c becomes in this case just 

c a x a2 ■ • • a n Q = [a^Q © (£>••• ® a n &\ ■ 
In the same way we have 

c bib*^ ■ ■ ■ b\n = [^n © © ■ ■ • © &;n] . 
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If we take now the inner product in the cyclic Fock space between these 
two vectors, then we get 

n-l 



p N (at ■ - ■ a n , 61 • • • b m ) = 5 nrn S ^{a l Vt, b* n+k Q) ■ ■ ■ (a n Q, b* 1+k Q) 

k=0 
n-l 

= 8 nm 2J y?(ai&„+fc) • • • <f(a n b 1+k ), 



k=0 

as required by our Definition 16.31 Thus the subalgebras A^, . . . , A^- 
are free with respect to (ip N ,p N ). The invariance assumption on the 
distribution with respect to (^Pn,Pn) is also easily verified and so we 
can apply our limit theorem. 

Let = <&(/) + ••• + <$(/)■ Since 

(cw(/i) • • ■w(/„)fi,cw(/ n+m ) • • ■uj(f n+1 )) cyc 

= Pat(5 , at(/i) • • ■ SN(fn), S'jvi/n+l) ' " - Sjv(/n+m)) 

we can take the limit as iV — > 00 and apply Theorem 16. 91 So it remains 
to identify the limits ip± and tp2 in the assumption of that theorem. One 
sees easily that 

Mfi,...,f n ) = lim N-<f N (q^(fi)---q%(f n )) 

(A,/ 2 ) ifn = 2 
otherwise 

and 

^(/l,--- ,f n ;9i,---,9m) 

= lim N ■ p N (qN(h) ■ tiffa) ■ ■ ■ ?ArOm)) 

(fi,9i) ifn = l = m 
otherwise 

This gives exactly our claim. □ 

7.2. Proof of Theorem 15. 71 From equation (J3J) we only have to prove 
that 

(40) (cp(di) • • -p(4)^, cp« +m ) • • •p« +1 )fi) cyc 

^ ^ 1pir(dl, ■ ■ ■ , <^n; • • • 5 ^n+m)- 
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Note that we can replace T> by D (g> L°°[0, 1], ip by ip ® r > where r is 
integration with respect to Lebesgue measure on [0,1], and for each 
iV G N, p(d) by 

p(d ® x(0, 1)) = pjv(rf) + z&(d) + ■ • • + p^(d), 
where we have put 

pjv(d) :=p(rf®x(4)) 
with x(In) denoting the characteristic function of the interval 

P fi - 1 i 



* \ N 1 N 

This fits into the framework of our general limit theorem by putting 
T — V, 

A N = A(V®L°° (0,1)), A l N =A{V®L co {I i N )), 

<p(a) = {afl,Q), p(a, b) = (c afl, c b*fl) cyc (a,&GA/v), 
and finally 

QM = P%(d) G 

One can check again by the same arguments as for the semi-circular 
case that A N , . . . ,A^ C An are free with respect to (<Pn,Pn)- Also 
the invariance assumption on the distribution with respect to (ipn, Pn) 
is easily verified. 

Since p(d) has, for each N, the same moments with respect to (fN 
and as Sjv(^) = QnW + ' ' ' + 5at(^)> we can calculate the moments 
of p(d) via SN{d) by sending N —>■ oo and invoking our limit theorem, 
Theorem 16.91 It only remains to identify the limits tpi and ip2 from the 
hypothesis of the theorem, and show that 

^2 VU^i, . . ■ , d n ; d 

a€NC(n,m) 

(41) = ^{d\,---,d n ,d 

n-\-l} • • • } dn-\-m) 

Tr£SNc(n,m) 

Note that each inner product appearing in the calculation of 

^N{p N {di)---p N {d n )) 

gives a factor 1/N; one inner product must be involved in any case 
to get a non- vanishing result, thus the sought limits single out exactly 
the contributions with one inner product. In the case of <^jv(Piv(^i) ' ' ' 
p l N (d n )) this means that p l N (d n ) must act as a creation operator, p N (di) 
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as an annihilation operator and all the other p's as preservation oper- 
ators, thus 

...,d n )= lim N ■ Lp N (p L N (d 1 ) ■ ■ -p l N (d n )) 

= (d 2 d 3 ■■•dn, d\) = i/j(did 2 ■ ■ ■ d n ) 

In the case of one has to note that the only relevant contributions 
to cp l N (di) ■ ■ •pj V (rf n )f2 are of the form: p l N (d n ) must act as creation 
operator; since c Q = 0, no annihilation operator can appear, but since 
c can also act by multiplication of arguments there might be a second 
action as creation operator (let's say of p l N (dk)), all the other p have 
to act as preservation operators. Thus the relevant contributions of 
cp l N (dx) ■ ■ •p % N {d n )Q, are the terms with k — 1, . . . , n of the form 

c(d l ---d k ® 4 + i ■ - -dn). 

Since we are looking for terms which give in the end exactly one inner 
product, the relevant action of c is given by multiplying arguments 
and yields terms of the form 

[dk+i ■ ■ ■ d n d\ ■ ■ - dk] for some k — 1, . . . , n. 

In the same way the relevant contributions of cp l N (d n+m ) ■ ■ ■ p l N (d n+1 )Q 
are of the form 

K+j-i • ■ ■ <+l<+rX+m-l • ■ • d *n+l) for some 1 = h ■ ■ ■ , ™.. 

Thus we have 

(42) tp 2 {di, . . ■ , d n ; d n +i, . . . 

= lim N ■ pN(p l N (di) ■ ■ -PN(d n ),p l N (d n+1 ) ■ ■ -p l N {d n+m )) 

N^oo 



= lim N ■ (cp^idi) ■ ■ -p l N (d n )n } cp z N (d* n+m ) ■ ■ -^« +1 )fi) cyc 

n m 

= 2J / J([^fc+i ' ' ' d n di - "dk], [d* n+ i_i ■ ■ ■ d* n+1 d* n+m ■ ■ ■ d* n+l ]) cyc 

k=l 1=1 

n m 

— ^(dn+i ' ' ' d n+m d n+ i ■ ■ ■ c? n +/-i4+i ' ' ' d n d]_ ■ ■ ■ dk) 

k=l 1=1 

Suppose a G NC(n,m) has more than one through-block. Then 

for each block B, tp B ( d i, ■ ■ ■ , d n ; d n+1 , d n+m ) = ^ B {du d n +m) 
by equation (|HUj). Thus 

ipa(di, • • • d n \ d n+ i, . . . j dn+m) — V , 7r(^l) • • • d n , <i n +l; • • • ) <^n+m) 

where 7r e Sj^c( n , m ) is the unique permutation whose cycle decompo- 
sition is the partition ir. 
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Figure 12. One of the 

possible cycles which 
contribute to the sum in ips i n 
the case B is a unique through 
block. 



Now suppose that a has only one through-block, let [a] be the set of 
all 7r G SNc( n , m ) whose cycle decomposition gives the partition o. If 
B is a block of a which is not a through-block then again by equation 
(|30|) ipB and ipB are equal. If B is the unique through-block then as in 
equation (p?T]) write B = {ji, . . . ,j r } U {j r +i, ■ ■ ■ , j r +s}- Then by (jUJ) 

ipB(dj x , ...,d jr ; d jr+1 , . . . , d jr+s ) = ip^d^, . . . , d jr ; d jr+1 , . . . , d jr+s ) 

— ^ ] ipc{dj 1 , . . . , c?j r , dj r+1 , . . . , dj r+s ) 

c 

where c runs over the cycles in 7r G [er] which give the block B. Hence 
iprj = Y^^eM V'tt and thus equation is proved. □ 

8. DlAGONALIZATION OF FLUCTUATIONS 

Let us now use our description of fluctuations of random matri- 
ces in terms of operators to diagonalize these fluctuations. The one- 
dimensional Gaussian case is well established in the physical and math- 
ematical literature (see, e.g., |Pol[ IA.IM| l.lohj ). whereas looking on 
the one-dimensional Wishart case and, in particular, on the multi- 
dimensional Gaussian case was initiated by Cabanal-Duvillard |G-D| . 
Indeed, trying to understand and reproduce the results of Cabanal- 
Duvillard was the original motivation for our investigations. 

Since the fluctuations are given by taking inner products in cyclic 
Fock space, we can achieve such a diagonalization by taking functions 
of our operators which yield elementary tensors in cyclic Fock space. 
This means we are looking for a kind of cyclic Wick products. 

8.1. Semi-circular case. We should look for cyclic analogues of the 
Wick products W(/i<8>- ■ -®f n )- Let us denote them by C(/i®- ■ -<g)/ n ). 
They should be determined by the property that 



cC(/i ® • • • <g> f n )n = [A ® • • • <g> /„]. 
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Notice that we have 
cW(f 1 ®---®f n )n = c(f 1 ®---®f n ) 

= [fl ® • • • <E> /„] + </i, fn) ■ C h ® • • • ® /n-1 

= c ® • • • (g) / n )ft + (/i, /„> • c W(/ 2 ® ... (g) 

thus we could define these cyclic Wick products by the following recur- 
sion: 

C{h ® ' • ' ® fn) = W{h ® • • • <g> / n ) - (A, /„) • W(/ 2 , . . . , /„_:). 

For n = 1, this means, of course, 

C(f) = W(f) = uj(f). 

If we put 

/:=/i = /2 = --- = /„ (with H/ll = 1), 
then we know that 

= w/)/2), 

where the {U n } n are the Chebyshev polynomials of the second kind. 
Let V n (x) = U n {x/2); then V n (u(f)) = W(f® n ). Now, if we write our 
cyclic Wick polynomials in this one-dimensional case as 

C(f*») = 2T n (u>(f)/2), 

then these T n must satisfy 

2T n = U n - U n - 2 in > 2) 

and 

T ± (x) = U 1 (x)/2 = x. 

This shows that the {T n } are Chebyshev polynomials of the first kind. 

Let us now consider the mult i- dimensional case. It is easy to see 
that if fi is orthogonal to f i+ i for alH = 1, . . . , k, then we have for all 
n(l), . . .,n(k) > that 

W{ff n{1) ®f® n{2) <g>. • -®/f n(fe) ) = W{f® n{l) )-W{f® n(2) ) ■ ■ ■ W{f® n{k) ). 

If we assume in addition that also f\ and f k are orthogonal then we 
get for the corresponding C: 

c(/f n(i) ® fT (2) ® • • • ®/r (fc) ) = n(i) ® / 2 ® n(2) ® • • • ® /r (fc) ) 

= n(1) ) • W(,/f n(2) ) • • • n(fc) ). 

The covariance between such functions in our random matrices is given 
by the inner product in the cyclic Fock space. If we have k, I > 2 and 
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fi, . . . , f k G Hm, and g u . . . , g x G Hr such that fa _L f i+1 for % = 1, . . . , A; 
and <7j _L for i = 1, . . . , I then we have 



lim k 2 Ti[V n{1) {X N (h)) ■ ■ ■ V n(k) (X N (f k )}, 

Tr[V m (i){X N (gi)) ■ ■ ■ V m Q)(X N (gi))) 

= (ur (l) ® • • • ® /r (fe) i, [<?r w ® • • • ® ^r (l) ])c yc . 

Thus we recover the results of Cabanal-Duvillard C-D for that case 



8.2. Compound Poisson case. Again, we are looking for polynomi- 
als C(di <g • • • <g d n ) which have the property 

cC(dt (g • • • <g d n )fi = [di <g • • • <g 

We have 

c <g • • • (g d n )Q = c (di (g • • • <g d n ) 

= [dt, ...,d n ] + [d n di, d 2 ,..., d n _i\ + ip(did n ) c (d 2 <g • • • <g d n _i) 
= cC{d\ (g • - • (g d n )fi + cC(d n di <g <i 2 <g ■ • • <g G? n _i)Q 
+ V(^i4)c W(<*2 (g • • • <g rf n _i)fi 

Thus we define the C's in the following recursive way: 

(43) W(di (g • • • <g d n ) = C(d 1 <g> • • • (g <Q 

+ C(d n rfi (g d 2 <g • • • (g d n -i) + ^(rfid„)^(rf 2 (g • • • (g 

There does not seem to be a nice closed form for this in the one- 
dimensional case. 

Let us also look at the mult i- dimensional situation. We model this 
by assuming that we have elements d\, . . . , d r G T> such that didj = 
for % ^ j. Then we have again for ^ i(j + 1) (j = 1, . . . , n) and 
fc(l),...,fc(ra) > that 

W{df^ ® . . . ® dgW) = • • • w{dfX ] )- 

If also z(l) 7^ then we have again equality between W and C, i.e. 

=w(^ 1) )--w(<5 b) ). 
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8.3. Poisson case. Let us specialize the general compound Poisson 
case to the usual Poisson case. 

The usual Poisson case is special within the class of compound ones 
by a very special state on T>. Restrict for the moment to one random 
matrix, i.e., the algebra D is generated by one element d. Then the 
fact that we have a free Poisson variable p(d) means that this d is a 
projection and thus 

*P{d k ) = ij){d) =: A. 

So we can identify p(d) = p(d 2 ) = . . . and everything reduces again 
to polynomials in just one variable p(d). Again one knows (see |Ansl 
Theorem 4.11]) that the linear Wick polynomials W n (d) := W(d® n ) are 
given by the orthogonal polynomials with respect to the distribution 
of p(d) (i.e. with respect to the Marchenko-Pastur = free Poisson 
distribution). Let us denote these polynomials by U n , then we have 

W n (d)=U n (p(d)). 

If we put C n (d) : = C(<i® n ), then the general relation between W and 
C becomes in this case: 

W n {d) = C n (d) + C n _ X (d) + \W n _ 2 . 

If we put C n (d) = T n (p(d)) for some polynomials T n , then the above 
tells us that 

n„ — xu n _ 2 = r„ + r n -i- 

This gives us exactly the polynomials {T n } which appear in Cabanal- 
Duvillard's results |C-D| . 

As an extension of this, we also get the multi-dimensional Poisson 
case: There the "diagonalizing polynomials" in more than one vari- 
able are given by alternating products in the one-dimensional linear 
polynomials {II n }. 

A more detailed investigation of this diagonalization of fluctuations 
will be presented in (KMS . 

9. Asymptotic freeness of Gaussian and constant 

matrices 

Our results about compound Wishart matrices can be considered 
as describing the limiting relation between Gaussian random matrices 
and constant matrices for special moments - namely those with pat- 
terns of the form X*D\XX*D2X ■ ■ ■ X*D n X. This raises, of course, 
the question whether we can say something substantial about the gen- 
eral relation between Gaussian and constant matrices. In view of the 
basic theorem of Voiculescu that Gaussian random matrices and con- 
stant matrices are asymptotically free, we would expect that we should 
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have the same kind of statement also on the level of fluctuations. We 
want to indicate here that this is indeed the case, thus providing strong 
evidence that our notion of "second order freeness" is indeed the cor- 
rect concept. Note that in the following definition we make a quite 
strong requirement on the vanishing of the higher order cumulants. 
This is however in accordance with the observation that in many cases 
the unnormalized traces converge to Gaussian random variables. Of 
course, if we have a non-probabilistic ensemble of constant matrices, 
then the only requirement is the convergence of ki; all other cumulants 
are automatically zero. 

Definition 9.1. 1) Let {A\, . . . , A s }n be a sequence of N x iV-random 
matrices. We say that they have a second order limit distribution if 
there exists a second order non-commutative probability space (A, <p, 
p) and G A such that for all polynomials p 1; p 2 , . . . in s non- 

commuting indeterminates we have 

(44) lim k l {tr[p l {A l ,...,A s )}} = l f(p 1 {a 1 ,...,a s )), 

(45) lim k 2 {Tr\ Pl (A 1 ,...,A a )],Tr\p 2 (A 1 ,...,A a )]} = 

p(pi(a 1 , . . . ,a s );p 2 (a 1 , . . .,a s )), 

and, for r > 3, 

(46) lim k r {Tr\p 1 (A l ,...,A s )],...,Tr\p r {A 1 ,...,A s )]} = 0. 

2) We say that two sequences of Nx iV-random matrices, {Ai, . . . , A S } N 
and {B 1 , . . . ,B t } N , are asymptotically free of second order if the se- 
quence {Ai, . . . , A s , Bi, . . . , B t } N has a second order limit distribution, 
given by (A, ip, p) and a±, . . . , a s , b±, . . . , bt G A, and if the unital alge- 
bras 

Ay := alg(l,a 1; . . . ,a s ) and A 2 := alg(l, &i, . . . , bt) 
are free with respect to (<f,p). 

Remark 9.2. Corollary 13.31 shows that a family {X N (f)}f €HR of Her- 
mitian Gaussian random matrices has a second order limit distribu- 
tion. Theorem 15 .31 identifies the limiting distribution in terms of cyclic 
Fock space, and in the proof of Theorem 15.31 we have in addition 
shown that the limiting distribution is free of second order in that 
if JCx, . . . , JC n C TC are orthogonal subspaces and A is the algebra gen- 
erated by {uj(f) | / G fCi} then Ai, . . . ,A n are free with respect to 
(ip,p) where ip(a) = (aQ,Q) and p(ai,a 2 ) = (caiQ,ca 2 Q) cyc . Thus 
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we have shown orthogonal families of Gaussian random matrices are 
asymptotically free of second order. 

Remark 9.3. Corollary 13.71 showed that if {X^}n is a sequence of 
complex Gaussian random matrices and P^{Di) = X^D^Xpf where 

{d[ N \ D2 , -D3 , • • • Dp }iv is a sequence of N x N complex matri- 
ces which converges in distribution to (d\, d 2) . . . d p ) in (T>, ip) then the 
family {P/v(_Dj)}j has a limiting distribution. Theorem 15.71 calculates 
the limiting distribution in terms of cyclic Fock space. In the proof 
of Theorem 15.71 we have shown that the limiting distribution is free of 
second order in that if didj = for i ^ j and Ai is the algebra gen- 
erated by p(di) then Ai, . . . ,A P are free with respect to (<p,p) where 
ip(a) = (afl,fl) and p(ai,a 2 ) = (c a%ft, c a^fl)^. Thus we have shown 
orthogonal families of Wishart random matrices are asymptotically free 
of second order. 

Now we can address the question of the relation between Gaussian 
random matrices and constant matrices. We can even be more general 
for the latter and consider random matrices which are independent 
from the Gaussian ones. 

Let, as usual, X N (f) (/ G TCm) be a family of Hermitian Gaussian 
random matrices 

X N (f) = (*„(/))£ r 

as in section 13.11 

Theorem 9.4. Let {Xjsr(f) \ f G Hm.}n be a sequence of Hermitian 
Gaussian N x N -random matrices and {Ai, . . . , A s }n a sequence of 
N x N -random matrices which has a second order limit distribution. If 
{Xjsr(f) I / G Hm}n and {Ai, . . . , A s }n are independent, then they are 
asymptotically free of second order. 

The proof of this theorem relies on the same kind of calculations as, 
for example, in |MNj . Since we do not want to go into random matrix 
calculations here, we defer more details about this to |KMSj . 

If the random matrices {Ai, . . . , A s } are non-random constant ma- 
trices with limiting distribution with respect to the trace, then all k r 
vanish identically for r > 2, thus they have a second order limit dis- 
tribution, and we get as a corollary of the above that the asymptotic 
freeness between Gaussian random matrices and constant matrices re- 
mains also true on the level of fluctuations, i.e., with respect to our 
concept of second order freeness. 
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A more systematic investigation of this concept will be pursued 
in forthcoming publications. In particular, fluctuations of Haar dis- 
tributed unitary random matrices from this point of view will be treated 
in PES]. 
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